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Abstract 

The main decoding algorithms for Reed-Solomon codes are based on 
a bivariate interpolation step, which is expensive in time complexity. Lot 
of interpolation methods were proposed in order to decrease the complex- 
ity of this procedure but stay expensive. Then Koetter, Ma and Vardy 
proposed in 2010 a technique, called re-encoding, which permits to reduce 
the running time. However this trick is devoted for the Koetter inter- 
polation algorithm. We propose a reformulation of the re-encoding for 
any interpolation methods. Finally, we apply it to the Welch-Berlekamp 
algorithm. 
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1 Introduction 

The algebraic decoding algorithms for the Reed-Solomon codes were very 
studied from the last decades, especially their decoding algorithms. The 
Welch-Berlekamp decoding method provides a simple approach to decode the 
Reed-Solomon codes up to the correction capacity of the code |WB86j . Then 
in 1997, Sudan generalizes this approach to decode beyond this bound, which 
supplies the first list decoding method for this family Suu97 . Two years 
latter, Guruswami and Sudan introduced another generalization of the last 
method to correct even more errors, that is up to the Johnson bound [GS99 . 
In these three previous methods, a bivariate interpolation step is needed, 
moreover their time complexities are given by this procedure which can be 
expensive. Thus a lot of algorithms were proposed to solve the bivariate inter- 
polation as efficient as possible [Koc96 ( GR06, AZ08, TrilO . Even with these 
computation improvements, the bivariate interpolation step is always expensive. 

In this way, Koetter, Ma, and Vardy introduced the notion of re-encoding 
|KMVll] , This trick does not decrease the asymptotic complexity in general, 
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but permits a considerable gain in practice. The re-encoding can be split 
into three phases as following : start to perform a translation by a codeword 
on the received word such that k positions become null, then modify the 
intern statement of the interpolation algorithm to have benefits of the null 
positions, finally remove after the interpolation the translation did at the first 
step. This technique implies to modify the intern state of the interpolation 
algorithm in relation with the null positions to speed up the running time of 
the interpolation step. This adjustement of the intern state is the main, and 
maybe the only one, drawback of the re-encoding. 

In this article, we propose a new reformulation of the re-encoding. This 
reformulation permits to use the re-encoding trick with any bivariate inter- 
polation algorithm without preliminary modification. However to be generic 
is under an assumption between the multiplicity and the K-degree of the 
interpolated polynomial. We apply this reformulation to the Welch-Berlekamp 
algorithm and we observe that the gain is huge. 

The organization of the paper is as follows. In the Section [2] we recall the 
main decoding algorithms based on interpolation as Welch-Berlekamp, Sudan 
and Guruswami-Sudan. The Section [3] is devoted to recall the principle of the 
original re-encoding and to introduce our reformulation. Finally, in the Section^] 
we apply our revisited re-encoding to Welsh-Berlekamp algorithm and present 
the performances. 

2 Interpolation based decoding algorithms 

2.1 Bivariate interpolation for the decoding 

Different decoding algorithms are based on the bivariate interpolation. This step 
is the most expensive and the asymptotic complexity is given by this bivariate 
interpolation. For example, Welch-Berlekamp, Sudan and Guruswami-Sudan 
algorithms are based on this procedure. Since the list decoding algorithm for 
alternant codes [ABCllj . is also based on interpolation step, we can also apply 
the re-encoding on it. In this paper we propose to deal only with the decoding 
algorithms for Reed-Solomon codes, it is why we propose firstly to recall the 
definition of this class of codes. 

Definition 1 (Reed-Solomon codes). Let ati, . . . , a n £ ¥ q be n distinct elements 
ofW q . The Reed-Solomon code of dimension k and support (a,) is given by 

RS(a, k) = {(P(ai), . . . , P(a„)) :Pe¥ q [X] <k }. 

The three following algorithms are based on the same principle: 

1. Compute a bivariate polynomial by interpolation of the received word y 
and the support of the Reed-Solomon code a. 
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2. Compute the univariate polynomial(s) P which generated the code- 
word^), as y-root(s) of the bivariate polynomial. 

The differences between the three following algorithm are the parameters of the 
bivariate interpolation and represent the most expensive cost in the complexity 
of these methods. In the following, we present quickly the main decoding algo- 
rithms for Reed-Solomon codes, the interested reader can find more information 
in [Barllj . 

2.1.1 Welch-Berlekamp 

The Welch-Berlekamp algorithm is an unambiguous decoding algorithm 
devoted to the Reed-Solomon codes WB86J. However, the most famous list 
decoding algorithms are based on this method. This is why, we propose to 
recall the main step of this algorithm. 

This method is based on the computation of the bivariate polynomial by 
interpolation satisfying 



(Wwb) 



O^Q(X,Y) = Qo(X) + YQ 1 (X), 
a J Q{a i ,y i ) = 0, Vi 6 {l,...,n}, 
degQo < n - 1 - 1, 
degQi < n — t — k, 



where t — [^V^J * s ^ ne correction capacity of the Reed-Solomon code. Thus we 
obtain the following pseudo-code: 

Algorithm 1: Welch-Berlekamp 
Input: The received word j/GFJ and the Reed-Solomon code C. 
Output: The codeword c £ C if it exists such that d(c, y) <t = [^^J , 
under the polynomial form. 

begin 

Q{X,Y) <- Interpolation(IP W s,C) 
return 



2.1.2 Sudan 

Sudan remarked that if we wish to correct more errors, with the Welch- 
Berlekamp algorithm, it could happen that there exist different Y"-roots of the 
bivariate polynomial satisfying the condition Sud97] . So he proposed to modify 
the interpolation problem in this way: 

(IPs)* < Q(a i ,y i ) = Q 1 Vte{l,...,n}, 

[ degQ^ <„_T-l-i(*-l), Vje{0,...,n- 
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We can summarize this method: 
Algorithm 2: Sudan 
Input: The received word y G F" and the Reed-Solomon code C. 
Output: A list of codewords Ci of C, such that Vi, d(v,Ci) < T. 

begin 

Q(X,Y) <r- Interpolation(IP s ,C) 
(Pi,...,Pl)<- Y-Roots(Q(X,Y)) 
Candidate <— {} 
for i G {1,...,£} do 

if d (Pi (a), y) < T then 
I Candidate -S— Candidate U {Pi (a)} 

_ return Candidate. 



2.1.3 Guruswami-Sudan 

Since the Guruswami-Sudan algorithm introduces the notion of root with multi- 
plicity from Sudan algorithm |GS99] . Let us to recall the definition of the Hasse 
derivative. 

Definition 2 (Hasse derivative). Let Q(X,Y) G ¥ q [X, Y] be a bivariate poly- 
nomial and a,b be two positive integers. The (a,b)-th Hasse derivative of Q 
is 

dc SxiQ) degy(Q) , . . 

Q [-,6] (x ,y)A £ (j{, l )'i-'- v ' :): '■ 

i=a j=b \ / V / 

Thanks to the Hasse derivative, we can give the definition of the root with 
multiplicity higher than one. 

Definition 3 (Root with multiplicity). Let Q(X,Y) £ ¥ q [X, Y) be a bivariate 
polynomial and (a,/3) € (F 9 ) 2 be a point. The point (a, (3) is a root with 
multiplicity s G N if and only if s is the largest integer such that for all i+j < s 

Q [l ' j] (a,f3) =0. 

Guruswami and Sudan noticed that it could happen that for some two poly- 
nomials Pi ,P 3o , we have y ko = P io (a k „) = P ja (a kQ ) and so the point (a ko ,y ko ) 
is a root of Q with multiplicity at least 2. So they proposed to add multiplicity 
constraint during the bivariate interpolation step. 

(IPgs) — < Q(a i: yi) = 0, with multiplicity s, Vi e {1, . . . , n} , 
I deg(Q,) < s(n-T) - 1 - j(k - 1), Vj G {0, ...,£}. 

thus the pseudo-code of this method is given by: 
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Algorithm 3: Guruswami-Sudan 



Input: The received word y G F™ and the Reed-Solomon code C. 
Output: A list of codewords Ci of C, such that Vi, d(v,Ci) < T. 

begin 

Q(X,Y) <- Interpolation(IP GS ,C) 
(P x ,...,P t ) <- Y-Roots(Q(X,Y)) 
Candidate <— {} 
for i G {1, . . . ,£} do 

if d (Pi (a), y) < T then 
|_ Candidate <— Candidate U {Pi (a)} 

_ return Candidate. 



2.2 Re-encoding 

Definition 4 (Interpolation problem). Let = {(ai, y%), . . . , (a„, ?/„)} C 
(F g x F q ) . The interpolation problem with multiplicity s associated to ^3, 
IP(?P, s), consists in finding Q(X,Y) such that the points (oti,yi) are a root 
of Q(X,Y) with multiplicity at least s. 

Lemma 1. Let s be an integer, a, (3 G F g and Q(X,Y) G F g [X, Y] a bivariate 
polynomial such that the point (a, /3) is a root of Q with multiplicity s. Then 
for all univariate polynomial P such that P{a) — ft, we have 

(X-a) s \Q(X,P(X)). 

Proof. See [Gur051 Lemma 6.6, p. 103]. □ 

We can generalize the previous lemma for all interpolation points, taking 
care the multiplicity. 

Proposition 1. Let^fi C (¥ q x¥ q ) n and s be a positive integer. The polynomial 
Q(X,Y) is a solution ofW($$,s) if and only if 

n 

V6 g {o, . . . , s - 1} , l[(x - ai y- b | qv>\x, l(x)), 

where Q^(X,Y) = Q[°' b l(X,r) is the b-th Hasse derivative in Y, and L(X) is 
the Lagrange polynomial of ^3, that is for all i G {1, . . . , n} , L(pti) — y{. 

Proof. See [AZ081 Proposition 1]. □ 

Let *P C (W q x V q ) n and Lk(X) be the Lagrange polynomial on k elements 
of *P, without lost in generality, assuming the k first positions. Let 



Then for alii G {1, . . . , k} , i\ = 0. 
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Proposition 2. Let C (F 9 x ¥ q ) n and s be a positive integer. The polynomial 
Q(X,Y) is a solution o/IP(Cp, s) if and only if Q(X,Y + Lk(X)) is a solution 
ofW($ n ,a). 

Proof. See |KMVlll Theorem 3]. □ 



3 Revisited re-encoding 

3.1 Re-encoding and interpolation algorithm 

A problem occurs with the re-encoding process: we have to modify the inter- 
polation algorithm in order to take care of the k first interpolation points in 
order to speed up the computation. So for each interpolation algorithm we have 
to adapt the initialization step to have the total benefits of the re-encoding 
step. As far we know, only the Koetter interpolation algorithm was modified 
to perform it. Although lot of interpolation algorithms were proposed, we can 
for the moment, use the re-encoding trick only using the Koetter interpolation 
algorithm. 



3.2 Revisited re-encoding 

Let L n (X) be the interpolation Lagrange polynomial of the set <p n . Thus \/i £ 
{1, . . . , n} , L n (cti) = Ti, and deg(L„) < n— 1. Since for all i £ {1, . . . , k) r{ = 0, 
it exists the polynomial L n -k(X) such that 

L n-k{X) = . -. 

Thanks to the previous remark on the Lagrange polynomials, we deduce the 
following proposition which is the key ingredient of our reformulation. 

Proposition 3. Let ^ n -k = {(otfc+i, L n - k (a k+1 )), •••,(«„, L n - k (a n ))} C 
(F, x ¥ q ) n - k , 

dcgy(fl) 

R(X,Y)= R i(X)Y j , 

be a bivariate polynomial over ¥ q and s be a positive integer such that s < 
deg F R. The polynomial R is a solution of IP(^3„_fc, s) if and only if 

degy(ii) / k \ 

Q(X,Y)= Yl [RjWHiX-aiy-nYi, 
is a solution o/lP(*p n ,s). 
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Proof. Since R is a solution of W(ty n -k, s), then for all b € {0, . . . , s — 1} 

n 

n (*-«*r b 

n 

n (^-«o s - b 

i=k+l 
n 

i=l 
n 

Since s > degy R,s — j > 0, the statement is hold. □ 

Thanks to the Proposition [21 we can compute a solution of the interpolation 
problem on <p = {(ati, yi), . . . , (a n , y n )} from <p„ = {(ai, 0), . . . , (a k , 0), 
(afe+i,?/fe + i -L fc (a fe+ i )),..., (a„,y„ - L fc (a rl ))}. Our revisited re-encoding 
could be seen as a decoding on the puncturing code. Since the Reed-Solomon 
code are MDS, the punctured code has the same dimension and is also a 
Reed-Solomon code. We could imagine to reiterate the re-encoding process 
taking *P„-fc = then the decoding will make on the multi puncturing code 
and the correction radius will increase. 

The Proposition [3] is under the assumption that the multiplicity s is greater 
or equal than the ^-degree of R, a solution of IP(Cp n _fc, s). Which is not a 
problem, because a solution of IP(*p„_fc,s + k), for all positive integer k, is 
also a solution of IP0P„_fc, s). However, this artificial augmentation of the 
multiplicity could increase also the X-degree of the solution, and so introduces 
some issue for the interpolation problem related to the decoding. This is why 
we deal only with the Welch-Berlekamp algorithm in Section |U 

4 Application to the Welch-Berlekamp algo- 
rithm 

4.1 Straightforward application 

In the Welch-Berlekamp decoding context, use the principle of the revisited re- 
encoding, is straightforward. Indeed, the only one condition in order to make 
practical our re-encoding is the multiplicity s is greater or equal than the Y- 
degree of the bivariate polynomial to compute. In the Welch-Berlekamp context 
the multiplicity s is exactly equal to the Y-degree, that is 1. Let S(X,Y) — 
S (X) + YSx(X) e ¥ q [X,Y] be a solution of IP(q3„_ fe , 1), then R given by the 



R^(X,L n ^ k (X)) 



E 

deg y (ii) 

E 



Ri(X)(L n - k {X)) 



3-b 



' k \ 

RjWjjiX-Oi)-* )(L n (x)y- b 
\ »=1 / 



Q(X,L n {X)). 
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Proposition [3] 



R(X,Y) = S (X)Y[(X-ai)+YSi(X), 

i=l 

is a solution of the IP(<p„, 1). Keeping the same notations and using the Propo- 
sition [5J we deduce directly a solution of the interpolation problem IP(^p, 1) 
from the simpler one IP(*p n _ fe , 1). Let Q(X, Y) G ¥ q [X, Y] such that 

Q{X,Y) = R(X,Y + L k (X)) 

= ^So(X)J[(X - at) + L k (X)Si(X)^ + YSx. 

In order to satisfy the interpolation conditions of the Welch-Berlekamp algo- 
rithm, we must have: deg S\ < n — t — k and deg So < n — t — 1 — k. It can be 
rewritten as 

Vj G {0, 1} , degS,- < n - t - k - 1 - j(-l). 

We deduce that the weighted-degree changes during the bivariate interpolation. 
Using the example describes below, we have to interpolate n — k points with the 
weighted-degree equal to -1, instead of to interpolate n points with the weighted- 
degree fe — 1, without modify the intern state of the interpolation algorithm. Let 
us to illustrate our claim by a toy example. 

Example 1. Let Fg, a be a 7-th primitive root of the unity such that 
a 3 + a + 1 = 0, C be the Reed-Solomon code RS((ct! l )i = o,...,6) 2) over 
Hence the Welch-Berlekamp method can correct up to I ^rpj = 2 errors. Let 
P(X) = a 6 X + a 5 G Fg[X] be the message under its polynomial form. The 
associated codeword is then (a, a 4 , a 6 , a 3 , a 2 , 1, 0). Assume there are 2 errors 
occur during the transmission in the first and 5-th positions, the received word 
is (a 5 , a 4 , a 6 , a 3 , a 3 , 1,0). 

Now let us to perform the revisited re-encoding. Using the previous notations, 
the Lagrange interpolation polynomial of the original interpolation points set 
<P„ = {(l,a 5 ),(a,a 4 ),(a 2 ,a 6 ),(a 3 ,a 3 ),(a 4 ,a 3 ),(a 5 ,l),(a 6 ,0)} is 

L n = X 6 + a 4 X 4 + a 2 X 3 + c?X 2 + a 2 X + a 2 . 

Assume that we want to vanish the 2 first points, then the Lagrange interpolation 
polynomial on these points is L k = a 4 X + 1, and the quotient 

L ^ = ,y W V = X " + a * x3 + X 2 + X + a. 

(X - 1){X - a) 

Then the new interpolation points set is 

¥ n _ k = {(a 2 , a 4 ), (a 3 , a 2 ), (a 4 , 0), (a 5 , a 6 ), (a 6 , a)} . 
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Hence the bivariate polynomial which interpolates ty n -k with multiplicity s = 1 
and weighted-degree -1 is 

S(X, Y) = Y( a 6 X 2 + a 4 X + a 3 ) + a 2 X + a 6 . 

=5i — Sq 

We deduce the polynomial which interpolates the = 
{(1, 0), (a, 0), (a 2 , a 4 ), (a 3 , a 2 ), (a 4 , 0), (a 5 , a 6 ), (a 6 , a)} , 

= YS 1 (X) + (X -l)(X -a)S Q (X) 

= Y(a 6 X 2 + a 4 X + a 3 ) + a 2 X 3 + aX 2 + a 5 X + 1. 

To finish the reconstruction step of the interpolation, we compute 

Q(X,Y) = R(X,Y + L n (X)) 

= Y(a 6 X 2 +a 4 X + a 3 ) + a 5 X 3 + a 6 X 2 + a. 

V v ' v ' 

=Qi =Qo 

In the Welch- Berlekamp algorithm the Y -root search is trivial. Indeed, it con- 
sists only in the division of the Qo by Qi 

P= -2a = a 6 X + a 5 , 
Qi 

which is exactly the sent message under the polynomial form. 
4.2 Performance 

Since it is one of the main goals of this article, we assume that we cannot modify 
the intern state of the interpolation algorithm. The asymptotic complexity of 
Koetter algorithm is 0(LN 2 ) where L is the F-degree of the bivariate polyno- 
mial Q and N the number of the linear constraints given by the interpolation 
conditions. Then the complexity of the standard Welch-Berlekamp algorithm is 
0{n 2 ). While the complexity of the original re-encoding is also 0(n 2 ), there is a 
non trivial speed up, since the first coordinates are 0, our revisited re-encoding 
exhibits an asymptotic complexity of 0((n— k) 2 ). We perform some timing tests 
with a naive implementation of the Welch-Berlekamp algorithm in the hight level 
computer algebra system: MAGMA [BCP97 . We propose to compare 3 decoding 
methods: Welch-Berlekamp algorithm without re-encoding, Welch-Berlekamp 
algorithm with the original re-encoding, and finally the Welch-Berlekamp with 
our revisited re-encoding. These 3 decoding methods were implemented with 
the same interpolation function, without modification or particular parameter- 
ization. The experimentations were done on a 2.13GHz Intel(R) Xeon(R). The 
timings presented in the Table [1] are in seconds unit for 100 iterations for each 
set of parameters. 
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TO 


C 


usual 


original re-encoding 


revisited re-encoding 


4 


RS[15,8] 


0.270 


0.230 


0.090 




RS[15, 10J 


0.260 


0.210 


0.080 




RS[15,12] 


0.250 


0.180 


0.050 




RS[15, 14] 


0.230 


0.180 


0.050 


5 


RS[31,16] 


0.930 


0.760 


0.250 




RS[31,20J 


0.820 


0.710 


0.220 




RS [31,24] 


0.820 


0.660 


0.100 




RS[31,28] 


0.840 


0.550 


0.080 


6 


RS[63,32] 


3.440 


3.130 


1.070 




RS[63,40J 


3.480 


2.890 


0.650 




RS[63,48] 


3.460 


2.580 


0.390 




RS[63,56J 


3.350 


2.300 


0.260 


7 


RS[127,64] 


16.760 


15.220 


4.440 




RS[127,80] 


16.840 


14.160 


2.570 




RS [127,96] 


17.400 


13.160 


1.260 




RS[127, 112] 


17.780 


11.600 


0.560 


8 


RS [255, 128] 


100.780 


92.070 


21.150 




RS [255, 160] 


104.100 


88.200 


11.300 




RS [255, 192] 


109.840 


83.910 


5.110 




RS [255, 224] 


113.550 


74.440 


1.950 



Table 1: Comparison between the Wclch-Bcrlekamp without re-encoding, with 
original re-encoding and revisited re-encoding. The shown timings are in second 
unit for 100 computations. 



5 Conclusion 

We introduce a new reformulation of the re-encoding process which permits to 
make it usable with any interpolation algorithm. However the assumption that 
the multiplicity s is smaller than the Y-degree is the price to be generic. We 
perform different tests with the Welch-Bcrlekamp algorithm showing that our 
reformulation provides an very important gain. 
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